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The 1-2 model 

Geoffrey R. Grimmett and Zhongyang Li 

Abstract. The current paper is a short review of rigorous results for the 1-2 
model. The 1-2 model on the hexagonal lattice is a model of statistical mechanics 
in which each vertex is constrained to have degree either 1 or 2. It was proposed in 
a study by Schwartz and Bruck of constrained coding systems, and is strongly con¬ 
nected to the dimer model on a decoration of the lattice, and to an enhanced Ising 
model and an associated polygon model on the graph derived from the hexagonal 
lattice by adding a further vertex in the middle of each edge. 

The general 1-2 model possesses three parameters a, b, c. The fundamental 
technique is to represent probabilities of interest as ratios of counts of dimer cover¬ 
ings of certain associated graphs, and to apply the Pfaffian method of Kasteleyn, 
Fisher, and Temperley. 

Of special interest is the existence (or not) of phase transitions. It turns out 
that all clusters of the infinite-volume limit are almost surely finite. On the other 
hand, the existence (with strictly positive probability) of infinite ‘homogeneous’ 
clusters, containing vertices of given type, depends on the values of the parameters. 

A further type of phase transition emerges in the study of the two-edge corre¬ 
lation function, and in this case the critical surface may be found explicitly. For 
instance, when a > b > c > 0, the surface given by ^/a = Vb + ^Jc is critical. 


1. Origin of the 1-2 model 

The 1-2 model originated in the work of computer scientists Schwartz and Bruck 
[23] on constrained coding systems. They studied an array of variables on the hexag¬ 
onal lattice El subject to the ‘not all equal’ constraint. Of particular interest to them 
was the asymptotic behaviour of the number of acceptable conhgurations on large 
bounded regions E[„ of the lattice, in the ‘thermodynamic limit’ as j' H. Using 
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the method of so-called ‘holographic reduction’, they were able to map their counting 
problem to one of counting the number of perfect matchings (or ‘dimer coverings’) 
on a certain graph derived from the hexagonal lattice. This last problem may be 
solved using the Pfaffian representation of Kasteleyn [15], Fisher [5], and Temperley 
and Fisher [24], 

When rephrased in the language of statistical mechanics, the work of Schwartz 
and Bruck amounts to the calculation of the partition function of the following 1- 
2 model of probability theory and mathematical physics. Let El = (V, E) be the 
hexagonal lattice of Figure 1.1, and let S = { — 1,1}® be the set of conhgurations 
of absent/present edges, where the local state —1 (respectively, 1) means absent 
(respectively, present). The sample space is the subset of E containing all cr G E 
such that: every vertex of El is incident to either one or two present edges. Thus, a 
conhguration comprises disjoint paths and cycles of present edges. 

Let EI„ be an n X n subgraph of El with periodic boundary conditions, and let 
Pn be the uniform probability measure on the set of 1-2 conhgurations on H^. We 
ask for properties of pn in the limit as n —)■ cxo. In particular, does the limit measure 
Poo := lim„^oo pn exist, and, if so, what can be said about the long-range correlations 
of edge-states under ^oo? It turns out that the connection to dimers may be exploited 
to answer such questions. 

The above system is a lattice model whose partition function can be computed 
by the calculation of certain determinants using the holographic algorithm of Valiant 
[25]. By introducing an invertible 2x2 matrix on edges of a graph and conducting 
a base change, the partition function of a general vertex model on a graph G is 
transformed into the partition function of perfect matchings on a certain decorated 
version of G. Valiant’s original holographic algorithm can be generalized by assigning 
different bases to different edges (see [16]), and the ensuing algorithm can be used 
to compute partition functions of a larger class of vertex models in polynomial time. 

The holographic algorithm seems, however, not to be the most efficient way to 
solve the 1-2 model. In particular, the correspondence between the 1-2 partition 
function and the dimer partition function on a corresponding Fisher graph, via the 
base change method, is not measure-preserving; thus, the computation of local statis¬ 
tics and related probabilities becomes complicated, even if possible. An alternative 
measure-preserving correspondence was introduced in [18], and this permits a num¬ 
ber of representations in closed form of probabilities associated with the 1-2 model. 
This method, and some of its consequences, will be described in the current review. 

Certain properties of the underlying hexagonal lattice are utilized heavily in this 
work, such as trivalence, planarity, and support of a 1? action. It may be possible 
to extend the results summarized here to certain other graphs with such properties, 
including the Archimedean (3,12^) and (4,8^) lattices. 
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Figure 1.1. The graph E[„ is an n x n lozenge wrapped onto a torus. 
A horizontal (respectively, NW/SE, NE/SW) edge is said to be of type 
a (respectively, type b, c). 


The formal dehnition of the 1-2 model is presented in Section 2. The model has 
strong connections to the dimer and Ising models as well as to a certain polygon 
model, and these connections are laid out in Section 3. Two approaches to the issue 
of phase transition are outlined in Section 4, using the geometry and the correlation 
structure, respectively, and an exact formula for the critical surface is given in the 
second case. 


2. Definition of the 1-2 model 

Whereas the 1-2 model of [23] is uniform in that there is only one parameter, we 
present here the more general three-parameter model of [16]. 

Let n > 1, and let ri, T 2 be the two shifts of El as in Figure 1.1. The pair 
(ri,r 2 ) generates a action on H, and we write = (14, for the (toroidal) 
quotient graph of El under the subgroup of Z^ generated by the powers r” and 
r^. The conhguration space is the set of all a G {—1,1}^" such that every 
n G is incident to either 1 or 2 edges e with cr(e) = 1. Note that a G if and 
only if —a G It is sometimes convenient to work with the vector a' given by 
a'{e) = i(l + a(e)). 

A vertex t G W is incident to three edges 61 , 62,63 of EI„, which are in the 
respective orientations: horizontal, NW/SE, and NE/SW. Let a G and let the 
signature at v be the triple (t|„ := (T'(e3)cr'(62)a‘'(ei) G {0,1}^ considered as a word 
with three letters in the alphabet with two letters. Let a, 6, c G [0, 00 ) be such that 
(a, 6, c) 4 (0,0,0). To the signature (t\^ is allocated the weight w{a\y) G {0,a,b,c} 
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000, 0 111, 0 001, a no, a 010, b 101, b 100, c Oil, c 


Figure 2.1. Each vertex has eight possible local conhgurations, six 
of which have nontrivial weights, a, b, or c as above. A vertex v is 
said to be of type s in conhgnration a if a\y has weight s G {a,b,c}. 
Whereas the type of an edge is deterministic (see Fignre 1.1), the type 
of a vertex is random. 

given in Fignre 2.1, and the weight fnnction in : —)■ [0, oo) is dehned by 

(2.1) w{a) = w{(T\y). 

Vn 


This gives rise to the probability measure /i„ given by 


(2,2) 


where 


(2.3) 

= Zn{a, b,c) = ^ w{a) 




is the partition function. A sample drawn (approximately) from pio is depicted in 
Figure 2.2. 

It turns out that the weak limit of the sequence (/i„ : n > 1) exists. However, no 
simple correlation inequality is known, and the proof of existence of the limit follows 
a different route using the relationship to dimer conhgurations outlined in Section 
3.1. 


Theorem 2.1 ([10, Thm 6.2]). The weak limit 

P '00 • lim Pyi 

n—>-oo 

exists and is translation invariant. 
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Figure 2.2. A realization of the (uniform) 1-2 model on Hio with 
a = b = c = 1, drawn by MCMC methods. Each vertex has degree 
either 1 or 2, and therefore all components are either paths or cycles. 


The inhnite-volume limit /ioo is of course a Gibbs state in the sense that it satishes 
the relevant Dobrushin-Lanford-Ruelle (DLR) condition (see [7, Sect. 4.4]). On the 
other hand, the structure of the space of such Gibbs measures is unknown. Neither 
is it known for which parameter-values //oo is ergodic (it is not ergodic under the 
conditions of [18, Thm 4.9] and Theorem 4.3(b) of the current paper). 

Remark 2.2. For edges e, f ofM. and for sufficiently large n, we write { 0^0 f)n 
for the two-edge correlation function of the measure /i„. By Theorem 2.1, the limit 

(aeCTf) = lim {aecrf)n 

n^oo 

exists. In Section f.3, we shall consider the limit hm|e_/|^.oo(c’‘eC’‘/) o,s an order 
parameter that is indicative of phase transition. 

3. The dimer, Ising, and polygon models 

The relationship between the 1-2 model and the dimer model is pivotal to the 
study of the former. Dimers are relevant also to the study of Ising models in two 
dimensions (see, for example, [17]), and the Ising model gives rise in turn to a ‘high 
temperature’ polygon model (see for example, [1, p. 75] and [9, 21, 26]). Therefore, 
the 1-2 model is connected hrmly to the Ising and polygon models. These connections 
play roles in the theory of the 1-2 model, and are summarized in this section. 
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Figure 2.3. Part of the hexagonal lattice El is drawn on the left, 
together with a 1-2 conhguration. The graph Ha on the right is ob¬ 
tained by replacing each face of H by a certain ‘gadget’. The left-hand 
1-2 conhguration gives rise to a dimer conhguration on the decorated 
graph on the right, as described in the text. The fundamental domain 
is outlined as a lozenge, and expanded in Figure 3.1. 

3.1. The dimer model. Let H = (V, E) be the hexagonal lattice, and let 
Ha = (Va, Fa) be the ‘decorated graph’ drawn on the right side of Figure 2.3. The 
graph in the hgure is obtained by replacing each face F of H by a certain ‘gadget’ 
comprising a path which is joined to the vertices of F in the manner drawn in the 
hgure. 

A 1-2 conhguration of the left side of Figure 2.3 gives rise to a dimer conhguration 
on the decorated graph on the right, as follows. Note hrst that V C Fa. Each v G V 
has three incident edges in Ha, and these edges in Ha may be regarded as the bisector 
edges of the three angles in H at v. 

Let cr G S. An edge e G Fa, incident to a vertex n G V, is designated present if 
and only if the two sides of the corresponding angle of H have the same states, that 
is, either both are present or both are absent. Once we have determined the states 
of the bisector edges of Ha, there is a unique extension to a dimer conhguration on 
Ha. See Figure 2.3. Note that the two 1-2 conhgurations a, —a G S give rise to the 
same dimer conhguration, and thus the above correspondence is two-to-one. 

Consider now the toroidal graph H„ and the corresponding decorated graph H„,a- 
The edges of H„ a are weighted, with the edge e = (i, j) having weight 

if e is a horizontal bisector edge, 
if e is a NW/SE bisector edge, 
if e is a NE/SW bisector edge, 
otherwise. 


(3.1) 


w. 




= 
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The weight of a dimer configuration is defined as the product of the weights of the 
edges that are present, and this gives rise (as in (2.2)) to a probability measure dn,A on 
dimer conhgurations. Now, 5„,A-probabihties may be represented as weighted counts 
of dimer configurations, and such quantities on planar graphs may be computed by 
the Pfaffian method of Kasteleyn [15], Fisher [5], and Temperley and Fisher [24], 
This leads to the following limit theorem. 

Theorem 3.1 ([18, Prop. 3.3]). Let a,b,c> 0. The limit measure 

5a := lim 5n,A 

n^oo 

exists and is translation-invariant and ergodic. 

Theorem 2.1 follows by Theorem 3.1 and the above correspondence between 1-2 
model conhgurations and dimer conhgurations. Note that the weak limit for 1-2 
measures need not be ergodic. 



Figure 3.1. A single fundamental domain of the decorated graph 
Eln.A obtained from the central lozenge of Figure 2.3. This illustrates 
the relationship between the fundamental domain and the original 
hexagonal lattice H. Note the homology cycles ly of the torus. 

A (Jn^A-probability may be expressed in terms of weighted counts of dimers, and 
these are studied via the Pfaffian representation. As explained in [10] and the refer¬ 
ences therein, the asymptotics (as n —)■ oo) of such Pfaffians depend on the so-called 
characteristic polynomial of the model. We do not dehne the characteristic polyno¬ 
mial here beyond saying that it is the determinant P(z, w) of the weighted adjacency 
matrix of the fundamental domain of Ha, oriented in a ‘clockwise odd’ manner, and 
illustrated in Figure 3.1. It is a function of the parameters a, b, c, and of two complex 
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variables z, w. It is shown in [18] that 

P{z, w) = + l)^ + + 6a^b^ + 6a^c^ + 6b^c^ — 2ab ^^ (a^ + b^ — c^) 

-2ac(w + -) (a^ + - b^) -2bc(- + -] (b^ + - a^) . 

\ w J ^ ^ \w z / ^ ^ 

The spectral curve is the zero locus of the characteristic polynomial, that is, the 
set of roots of the equation P{z,w) = 0. As explained in [18], it is important to 
understand the intersection of the spectral curve with the unit torus 

= {{z,w) e : l^l = Itcl = l}. 

It turns out that the intersection is either empty or is a single real point (1,1). 
Moreover, when P(l, 1) = 0, the zero (1,1) has multiplicity 2. Evidently, 

(3.2) P(l, 1) = (a^ + - 2ah - 2bc - 2ac)^ 

and therefore the spectral curve intersects if and only if 

(3.3) \/a ± \/b ± \/c = 0. 

We shall return to this equation in the study of phase transition in Section 4.3. 


3.2. The half-edge graph. When considering correlations, it will be conve¬ 
nient to work on a graph derived from the hexagonal lattice by replacing each edge 
by two half-edges. Let AEI^ = (AVn, AEn) be the graph derived from EI^ = (14, En) 
by adding a vertex at the midpoint of each edge in En- Let MEn = {Me : e G be 
the set of such midpoints, and AVn = 14 U MEn- The edges AEn are the half-edges 
of En, each being of the form {v, Me) for some v E Vn and incident edge e E En- 
The 1-2 model on EI^ can be viewed as a spin-model on the set MEn of midpoints, 
as we explain next. 


3.3. The Ising model. It turns out that the 1-2 model is the marginal of a 
certain Ising-type model on the half-edge graph AEI„ of Section 3.2, that is reminis¬ 
cent of the Edwards-Sokal coupling of the Potts and random-cluster measures (see 
[7, Sect. 1.4]). It is constructed via a weight function on conhguration space, using 
weights that are permitted in general to be C-valued. 

Let E® = {—1, and 'E' = {—1, -|-1}^". An edge e G Pn is identihed with 

the element of MEn at its centre. A spin-vector is a pair (cr®, a'') G E® x S'" with 
(T® = {a^^s '- V E Vn, s = a, b, c) and a'' = {ay : v G 14), where ay^a, (^y,b, (^v,c denote 
the spins on midpoints of the corresponding edges of given types incident to n G 14 
(see Figure 1.1). We allocate the (possibly complex) weights 

(3.4) w{a'', cr®) = n (1 + ea<^vO-y,a){^ + e6cr^cr„,b)(l + ecO‘^(L„,c), 

VGVn 
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where ea,eb,ec G C are constants associated with horizontal, NW/SE, and NE/SW 
edges, respectively. In (3.4), each factor 1 -|- ega^ay^s (s = a,b,c) corresponds to a 
half-edge of H^. 

When considering the relationship to the 1-2 model, it will be convenient to 
choose the constants e*, s = a, b, c, as follows. Let a,b,c> 0, and 


(3.5) 


A 


a — b — c 
a + b + c' 


B 


b — a — c ^ c — a — b 
a + b + c' a + b + c' 


where we assume for simplicity that ABC ^ 0. The appropriate values of the are 


(3.6) 



Theorem 3.2 ([10, Sect. 4]). 

(a) Marginal on 14. For a'' G S'", w{a'') := satisfies 

(3.7) = 2l®"l JJ (l -F elau(Ty) . 

g={u,v)&En 

That is, the marginal weights on YA are those of an Ising-type model on EI^ 
with (possibly complex) edge-interactions. Here, Cg denotes the parameter 
associated with edge g = {u,v). 

(b) Marginal on MEn- Let a,b,c > 0 and assume the satisfy (3.5)-(3.6) 

where ABC 4 0. For cr® G S®, w{a^) := ^ satisfies 

R'(<T®) =11(1 + Aa^^bCr^^c + 

veVn 

OC /in(cr4- 

That is, the marginal weights on S® are proportional to those of the 1-2 
model. 

(c) Two-edge correlation function. Let e = {u,v),f = {x,y) G Subject to 
the notation of part (b) above, the two-edge correlation function of Remark 
2.2 satisfies 


(3,8) 


(ffeff/)., = E De,f(<T'')w(a'') / Yi 


w[a 


xVgS'' 


T-VgS'' 


where 


p, / v\ ^e(^u + + (^y) 

) = 


{l + el){l + e}) 
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If the weights w{a'^) of (3.7) are real and non-negative (which they are not in 
general), the ratio on the right side of (3.8) may be interpreted as an expectation. 
The weights w{a'^) correspond to a ferromagnetic Ising model if and only the edge- 
weights of (3.7) satisfy 0 < < 1. Unfortunately, this never occurs with the eg 

derived from the 1-2 model as in (3.5)-(3.6). If, however, one assumes that the a, b, 
c satisfy the acute angle condition 

2, 1,2 I 2 72 ,2, 2 2, 2, 7,2 

a < b + c , b < c + a , c < a + b , 

then —1 < Cg < 0, and the corresponding Ising model is antiferromagnetic. Since 
HI is bipartite, this process may be transformed into a ferromagnetic system by 
changing the sign of every other vertex (see [ 6 , p. 17]), and this transformation 
greatly facilitates its analysis. 

The above Ising model may be regarded as a special case of the eight-vertex 
model of Lin and Wu [20]. 

3.4. The hexagonal polygon model. Let ]H[„ = (Vn, En) as before, and let 
Yin = {0,1}^". The sample space of the polygon model on ]H[„ is the subset 11 ])°^ c 
n„ containing all tt = ( 7 r(e) : e G En) G Tin such that 

(3.9) ^ 7r(e) is either 0 or 2, n G I 4 . 

eBv 

Each TT G nP°h j^iay be considered as a union of vertex-disjoint cycles of ]H[„, together 
with isolated vertices. We identify tt G with the set {e G : 7 r(e) = 1} of ‘open’ 
edges under tt. Thus (3.9) requires that every vertex is incident to an even number 
of open edges. 

Let Ca, Cfe, Cc G C. To the conhguration tt G 11 ])°^^ assign the (possibly complex) 
weight 

(3.10) W{7r) = e2k(a)|^2|7r(6)|^2K(c)|^ 

where 7 r(s) is the set of open s-type edges. The weight function w gives rise to the 
partition function 

Let a,b,c > 0. We now choose the constants Cs, s = a, b, c, according to (3.5)- 
(3.6), where it is assumed that ABC 7 ^ 0. The corresponding polygon model is 
related to the high-temperature expansion of the Ising-type model of Section 3.3 
(see, for example, [ 8 , eqn 5.1] and [9, Thm 1.7]). 

In considering correlation functions, it is convenient to view the polygon model 
on the half-edge graph AM.n of Section 3.2. A polygon conhguration on ]H[„ induces a 
polygon conhguration on A]H[„, namely a subset of AE^ with the property that every 
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vertex in AVn has even degree. For an a-type edge e G the two half-edges of e 
have weight each (and similarly for b- and c-type edges). The weight function w 
of (3.10) may now be expressed as 


w{7r) 


-|7r(a)| k(fe)| |7r(c)| 
a '-f) '-C 1 


where 7r(s) is the set of open half-edges of type s, as tt ranges over polygon conhgu- 
rations on A]HI„. 

Let e, / G MEn be distinct midpoints of and let Ilej be the subset of all 

TT G {0,1}^®" such that: (i) every v G AVn with v ^ e, f is incident to an even 
number of open half-edges, and (ii) the midpoints e and / are incident to exactly 
one open half-edge. Let 

(3-11) V„(e,/) = %gi, 

where 

(3.12) Zn,e^f.= Y. 

TTGlle,/ 


Theorem 3.3 ([11]). Subject to (3.6), the two-edge correlation function {ae<7f)n 
of the 1-2 model on EI^ satisfies {ae(Tf)n = Mn{e, /). 

The polygon model with general parameters is studied in [11]. 


4. Phase transition 

We discuss two forms of phase transition for the 1-2 model. Of these, the hrst 
concerns the existence (or not) of inhnite ‘homogeneous’ clusters of El containing 
vertices of the same type, and the second considers as order parameter the limiting 
two-edge correlation function. Thus, the hrst studies the geometry of the model, and 
the second its correlation structure. There may exist other forms of phase transition, 
as yet unstudied. 

4.1. Occurrence of paths. Every connected component (‘cluster’) in a real¬ 
ization of the 1-2 model is either a self-avoiding path or a cycle (see Figure 2.2). It 
turns out that all such clusters are /ioo-a-s. hnite, when a,b,c> 0. We formalize this 
statement in this subsection, and begin with an exact formula. Using the correspon¬ 
dence between 1-2 model conhgurations on El and dimer conhgurations on Ha, as 
described in Section 3.1, we have the following. 

Theorem 4.1 ([18, Thm 3.4]). Let a,b,c> 0 and let //oo be the limit 1-2 measure 
of Theorem 2.1. Let i be a self-avoiding path o/H„ containing I -\- 1 edges, and write 
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Ei = {ck = {uk,Vk) : 1 < k < 1} for the set of bisector edges o/H a encountered 
along 1, as in Figure 2.3. Then 


/ioo(^ is present) = ^ j |Pf ^1 > 

where We is the weight of the edge e in Ha, Kff^ is the submatrix of the inverse of the 
weighted adjacency matrix o/Ha with rows and columns indexed by ui,vi,... ,ui,vi, 
and Pf M is the Pfaffian of the matrix M . 

The mass transport principle, introduced in [3, 13], is a valuable tool in the 
study of interacting systems including percolation and self-avoiding walks on Cayley 
graphs, see [2, 12, 14], It may also be applied to the 1-2 model, where it is used to 
prove the following. 

Theorem 4.2 ([18, Thm 2.4]). Ifa,b,c > 0, we have that 
^oo (there exists an inhnite path) = 0. 

4.2. Existence of infinite homogeneous clusters. The concept of ‘phase 
transition’ hinges on the non-smoothness of some so-called ‘order parameter’. For 
the Ising model, one may take as order parameter the magnetization at the origin 
in the inhnite-volume measure with -|- boundary conditions. This corresponds in the 
universe of percolation and the random-cluster model (see [7]) to studying whether 
or not there there exists an inhnite cluster. By Theorem 4.2, the 1-2 model possesses 
no inhnite cluster for any values of a,b,c > 0. ‘Clusters’ may however be dehned in 
another manner. 

Let a E Ti. Each vertex of H has a random type, given in Figure 2.1. For 
s G {a, b,c}, a types cluster is a maximal connected subgraph of H every vertex of 
which has type s. As illustrated in the hgure, type s comes in two sub-types; for 
example, a type-a vertex has signature either 001 or 110. Thus, one may speak of 
a 001-cluster, etc. By examining the hgure again, it is seen that a type-a cluster 
is either a 001-cluster or a 110-cluster, but may not contain vertices of both types 
(and similarly for types b and c). A homogeneous cluster of cr G S is a tc-cluster for 
some w G {0,1}^, w ^ 000, 111. We concentrate now on the existence (or not) of an 
inhnite homogeneous cluster. 

Theorem 4.3. Let a,b,c> 0. 

(a) ([19, Thm 1.1]) Let w G {0,1}^, w ^ 000,111. The number of infinite 
w-clusters is /ioo-o-S- no greater than 1. 
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(b) ([18, Thm 4.4, Prop. 4.7]) Fix 6, c > 0. For sufficiently small a > 0, there 
exists poo-o-S- no infinite type-a cluster. For sufficiently large a, the poo- 
probability that the origin belongs to an infinite type-a cluster is strictly 
positive. 

Part (a) is proved in [19] using an adaptation of the method of Burton and Keane 
(see [4, 22]) to the 1-2 model, subject to the complication that poo does not have 
the so-called ‘hnite energy property’. It is unknown whether inhnite tc-clusters and 
tc'-clusters can coexist with w ^ w'. 

Part (b) indicates the existence of a phase transition. It is not known if there 
exists a single critical point Oc = ac{b,c) for the given property. Furthermore, there 
is currently no indication of the exact value of such a point. 

4.3. Non-analyticity of the two-edge correlation function. For distinct 
edges e, / G E, we write 

(cTeCr/) = lim {aeaf)n, 

n—>-oo 

which exists by Theorem 2.1, (see [10, Thm 6.2]). We consider here the asymptotic 
behaviour of ffieCTf) as \e — f \ —?■ oo. The behaviour of this limit is unknown in 
general, but a great deal is known if e and / are related in the ‘diagonal’ manner of 
the forthcoming assumption (4.1), as illustrated in Figure 4.1. 

Theorem 4.4 ([10, Thm 3.1]). Let a,b,c> 0, and let e, / G E 6e NW/SE edges 
such that: 

there exists a path i = i{e, f) o/H^ from e to f 

(4.1) 

using only horizontal and NW/SE half-edges. 

(a) Let a> b > 0. For c > 0 satisfying 

either \/a > Vb + \/c or \/c > \/a -|- Vb 

except possibly on some set of isolated points, the limit hm|e_/|^oo(c’‘eC’‘/)^ 
exists and is non-zero. 

(b) If a > b > 0 and 

ffia — Vb < \/c < y/a + Vb, 
then ffieCTf) 0 as \e — f\ ^ oo. 

By Theorem 4.4, when a > 6 > c > 0, the phase transition occurs when ffia = 
y/b + \fc. The proof is along the following lines. The square of the two- 

edge correlation may be expressed as the determinant of an explicit block Toeplitz 
matrix (this is where (4.1) is used); see the forthcoming Lemma 4.5. Its limit as 
\e — f\ —)■ oo is given by Widom’s theorem (see [27, 28] and [10, Thm 8.7]) as 
the determinant of the limiting (inhnite) block Toeplitz matrix. This determinant 
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Figure 4.1. A path i comprising horizontal and NW/SE mid-edges, 
connecting the midpoints of two NW/SE edges e and /. 


is complex analytic with respect to the parameters a, b, c except when the spectral 
curve has a unique real zero on the unit torus. As remarked at the end of Section 
3.1, the last occurs if and only if ^/a ± Vb ± ^/c = 0. When a > b > c, this equation 
becomes -s/a — y/b — y/c = 0. 

The ‘isolated points’ of part (b) arise through the use in the proof of the fact 
that, for an analytic function A, either A = 0 or the zeros of A are isolated. 

It is easily seen from Figures 2.1 and 4.1 that 


{(TeCTf) = 1 if either a,b > 0, c = 0, or a = fo = 0, c > 0, 


and part (b) of Theorem 4.4 follows by the analyticity. For part (a), one uses the 
representation of the 1-2 model as the Ising model of Section 3.3. 

The key step in the proof of Theorem 4.4 is the following exact formula. Let Yi 
be the 2x2 matrix 


Fi = 






and let Y 2 k be the 4/c x Ak block diagonal matrix with diagonal 2x2 blocks equal to 
Yi. That is. 




0 •• 

■ 


0 

Fi •• 

■ 0 

Y2k = 




Vo 

0 • 

■ yJ 


Lemma 4.5 ([10, Lemma 8.4]). Suppose the path i = i{e,f) of (4.1) passes 2k 
bisector edges o/Ha. We have that 


(UeU/) =Pf [Y^k + ^cK^^ 


where ^ is as in Theorem f.l. 












THE 1-2 MODEL 
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